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We study the elementary excitations of Bose-Einstein condensates in a one- 
dimensional periodic potential and discuss the stability of superfluid flow 
based on the Kronig-Penney model. We analytically solve the Bogoliubov 
equations and calculate the excitation spectrum. The Landau and dynamical 
instabilities occur in the first condensate band when the superfluid velocity 
exceeds certain critical values, which agrees with the result of condensates in 
a sinusoidal potential. It is found that the onset of the Landau instability 
coincides with the point where the perfect transmission of low- energy excita- 
tions is forbidden, while the dynamical instability occurs when the effective 
mass is negative. It is well known that the condensate band has a peculiar 
structure called swallowtail when the periodic potential is shallow compared 
to the mean field energy. We find that the upper side of the swallowtail 
is dynamically unstable although the excitations have the linear dispersion 
reflecting the positive effective mass. 

PACS numbers: 03.75. Lm, 05.30. Jp. 03.75. Kk 



1. INTRODUCTION 

The realization of Bose-Einstein condensates (BECs) in a optical lattice 
has attracted much attention of both theorists and experimentalists ^. One 
of the interesting problems in this system is the breakdown of the superflu- 
idity. In a recent experiment, the stability of BECs in a moving ID optical 
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lattice has been investigated, it has been demonstrated that the Landau and 
dynamical instabilities destabilize the system when the velocity of the lat- 
tice exceeds certain critical values^. In addition, in a number of theoretical 
works, the stability of BECs in an optical lattice have been studied in great 
detail EHEEEEI 

In the present paper, by means of the Kronig-Penney (KP) model, we 
study the elementary excitations and the stability of BECs with supercurrent 
in a ID periodic potential. In our previous work, the excitation spectrum 
of current-free BECs in a Kronig-Penney potential has been calculated by 
applying the Bloch theorem to the solution of Bogoliubov equations for a 
single-barrier problem^. We will extend this theory to the case of a current- 
carrying BECs. 



2. MEAN FIELD THEORY AND CONDENSATE BAND 

We consider a quasi-lD BEC at the absolute zero of temperature in a 
Kronig-Penney potential, which is a periodic array of 5-function potential 
barriers, 



V(x)=V J2 S(x-ja) 



(1) 



where a is the lattice constant, and Vq is the potential strength. We consider 
that the condensate has supercurrent flowing through the periodic potential, 
which corresponds to the situation where the optical lattice created by two 
counter-propagating laser beams with frequency difference is moving at a 
constant velocity^. 

Our formulation of the problem is based on the mean field theory, which 
consists of the time-independent Gross-Pitaevskii equation and the Bogoli- 
ubov equations Assuming that the condensate is weakly perturbed and 
seeking the solution of the time-dependent Gross-Pitaevskii equation in the 
form ty(x,t) = ^o( x ) + u(x)e~~ — v(x)*e~ , one obtains them, 



2m dx 2 



and 



*2 



H + V(x)+g\V (x) 



-H 



*o(x) = 0, 




2m dx 2 



li + V{x)+2g\* Q \ 



(2) 

(3) 
(4) 
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where is the chemical potential and g is the coupling constant of the 
interparticle interaction. The excitation energy e tells us the stability of a 
BEC. The appearance of the negative excitation energy is a signature of the 
Landau instability, while the appearance of the complex excitation energy 
is that of the dynamical instability. 

At first, we shall solve Eq. ©. We assume that the condensate stands 
in the first condensate band. Then, assuming that the lattice constant a is 
much larger than the healing length £, one can approximately obtain the 
solution of Eq. (J2J) as 



~ g 3eQ (tanh ( |X -^ I+X °) - sgn(x)i^) e^-^o), 

' ^ a < x < ( j + a, (5) 



*o(z) 



2 J V 2 

where en = Tfcr and £ = , Constants xq and 9o are determined by 

the boundary conditions at x = ja. 

Imposing the Bloch theorem ^o(x + a) = ^o(x)e~ on the condensate 
wave function ®, we can calculate the first condensate band E, the group 

velocity v g = Jj| and the effective mass m* = fj^J as functions of the 
condensate quasimomentum K as shown in Fig. The swallowtail structure 
in the first condensate band exists around the edge of the first Brillouin zone 
(Fig. ID^a)), which disappears when the potential strength is sufficiently large 
compared to the mean field energy^. Several authors have pointed out that 
the upper side of the swallowtail is dynamically unstable* 4 * 1 ^ 14 ! On the other 
hand, the system is expected to be stable for the long wavelength phonon 
in the upper side, because the effective mass is positive (Fig. m c ))- We will 
confirm these two fact, namely the dynamical instability and the existence 
of the stable phonons in the upper side. 




Fig. 1. First condensate band (a), group velocity (b), and effective mass (c) 
with (a, V ) = (15£ av , 5sm av £ av ), where £ av = j== , c s = -\/^fS n av = ^p, 
and Nq is the number of condensate atoms. 
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3. EXCITATION SPECTRUM AND STABILITY 



Next, we analytically calculate the excitation spectrum. We shall first 
find solutions of the Bogoliubov equations © in the region \x\ < # where 
only a single barrier exists. There exists two independent solutions with the 
excitation energy e of the single barrier problem, corresponding to two types 
of scattering process. One is the process in which an excitation comes from 
the left , and the other from the right (ip r (x)). ip l (x) is written as^ 



^j\x) 



U1 )+r 1 r 2 wr 3 



\ Vl 



+ 6 1 



U4 



x < 0, 



x > 0, 



(6) 



where the coefficients r , a 1 , r, and 6 1 are the amplitudes of the reflected, the 
left localized, the transmitted, and the right localized components, respec- 
tively. All the coefficients can be determined by the boundary conditions at 
x = 0. Once the coefficients of ip^(x) is obtained, the coefficients r r , a r , t r , 
and b r of ip T (x) can be easily obtained as well^M The transmission coefficient 
|t'| 2 strongly depends on the supercurrent as follows. If there is no supercur- 
rent, the transmission coefficient has a kind of resonant behavior with a peak 
at e = 0, which is called anomalous tunneling^^. As the supercurrent in- 
creases, the width of the peak becomes narrow; consequently, the anomalous 
tunneling disappears at a certain critical value of the supercurrent^. 

We can write a general solution of the Bogoliubov equations (j3J) in 
| x | *C 7y as a linear combination of ip l (x) and ip T (x): 



ip(x) = aij/(x) + 0ip T (x), 



i i a 
M < 2> 



(7) 



where a and (3 are arbitrary constants. Extending this solution to all regions 
of x by means of the Bloch theorem, one can obtain the relation between 
the energy e and the quasimomentum q of the excitation: 



cxp 



(-k\ + k 2 )a 



t exp 



2h 

.(-2q + k 1 + k 2 )a 



2Ti 



+ (t l f - rV r )exp 
+ t r exp 



1 2^ . 
. (2q -ki- k 2 )a 
1 2h 



(8) 



where ki }2 satisfies the Bogoliubov excitation spectrum in a uniform sys- 
tem^, 
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Fig. 2. (a) Stability phase diagram of the condensate in a KP potential, 
where (a, Vq) = (15£ av , 5<?n av £av)- The light shaded area and the dark shaded 
area correspond to the regions of the Landau instability and the dynamical 
instability, respectively. This figure is focused to a region close to the swal- 
lowtail, (b) Second band of the excitation spectrum. 



As well as the case of a current-free condensate, we clearly see from Eq. (JSJ) 
that the tunneling properties in the single barrier problem determine the 
band structure of the excitation spectrum. 

Now let us discuss the stability of condensates in a KP potential. We 
note in advance that there are three important values of the condensate 
quasimomentum, Km and Ke- Ka, corresponds to the point where the 
anomalous tunneling disappears, Km to the point where the group velocity 
reaches its maximum value and the sign of the effective mass changes, and 
Ke to the edge of the swallowtail, respectively. 

Calculating the excitation spectrum by solving Eq. Q, we obtain the 
stability phase diagram as shown in Fig. |2fa). The light (dark) shaded area 
corresponds to the region of the Landau (dynamical) instability where the 
negative-energy (complex-energy) excitation exists. Both of the instabilities 
appear in the first band of the excitation spectrum. In the lower side of the 
swallowtail, K = Ka is the onset of the Landau instability around q = and 
the dynamical instability around q = ^ , while K = Km is the onset of the 
dynamical instability around q = 0. In the upper side, the system is always 
dynamically unstable, which agree with the result of several recent works. 
On the other hand, as shown in Fig. Efb), there also exists the gapless and 
linear dispersion in the upper portion, reflecting the positive effective mass. 
As K in the lower portion of the swallowtail approaches Ke, the bottom of 
the second band of the excitation spectrum approaches the origin. When K 
reaches Ke, the second band turns into a gapless dispersion. Consequently, 
there exists the gapless and linear dispersion in the upper portion, reflecting 
the positive effective mass, even though the system is dynamically unstable 
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due to the first Bogoliubov band. 



4. CONCLUSION 

In summary, we have investigated the excitation spectrum and the sta- 
bility of BECs in a KP potential. We have found that the onset of the 
Landau instability coincides with the point where the anomalous tunneling 
disappears, while all the excitations in the first band of the excitation spec- 
trum exhibits the dynamical instability when the effective mass is negative. 
It has been shown that the second band of the excitation spectrum takes 
the gapless and linear form for the long-wavelength excitations in the upper 
side of the swallowtail, while the system is always dynamically unstable. 
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